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Abstract 

An adaptive unstructured meshing technique for 
mechanical and thermal stress analysis of built-up 
structures has been developed. A triangular membrane 
finite element and a new plate bending element are 
evaluated on a panel with a circular cutout and a frame 
stiffened panel. The adaptive unstructured meshing 
technique, without a priori knowledge of the solution to 
the problem, generates clustered elements only where 
needed. An Improved solution accuracy Is obtained at a 
reduced problem size and analysis computational time 
as compared to the results produced by the standard 
finite element procedure. 


Nomenclature 

finite element area 
strain-displacement matrix 
material stiffness matrix 
bending rigidity 
modulus of elasticity 
load vector 
nodal spacing 
stiffness matrix 
bending moment 

element Interpolation function matrix 
pressure 

vector of stress components 
temperature 

reference temperature for zero stress 
element thickness 
displacement components 
local coordinate directions 
global coordinate directions 
coefficient of thermal expansion 
vector of coefficients of thermal expansion 
vector of element nodal unknowns 
vector of strain components 
eigenvalue 
Poisson's ratio 
key parameter for remeshlng 
stress components 
principal stress components 

‘Aerospace Technologist, Aerothermal Loads Branch, Structural 
Mechanics Division 
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Subscripts 
b bending 

m membrane 

p pressure 

T thermal 

Superscripts 
T transpose 


Introduction 

Accurate prediction of the structural response 
induced by both mechanical and thermal loads Is an 
important factor In the design of many structures. Intense 
pressure and heat transfer rates may produce severe 
stresses that reduce the structural performance and may 
cause structural failure. Research is underway to 
Improve the efficiency and accuracy of the structural 
analysis procedure. 

An adaptive unstructured meshing technique 1 has 
been shown to Improve the efficiency and accuracy^ of 
the high speed How analysis, the thermal, and structural 
analyses of continuum structures In two dimensions by 
the finite element method. The technique generates a 
new mesh based on the solution obtained from an earlier 
mesh. The new mesh consists of small elements in the 
regions with large changes in solution gradients and 
large elements In other regions where the gradient 
changes are small. Because the technique generates 
proper element sizes automatically, it Is especially 
suitable for complex problems where a priori knowledge 
of the solutions does not exist. 

Herein, the technique has been extended to 
construct adaptive meshes for structural analysis of ‘built- 
up* structures. Such structures are commonly modelled 
by using two-dimensional membrane and plate bending 
finite elements. A new triangular plate bending element 
to be used with the adaptive unstructured meshes for 
efficient thermal stress solution is Introduced In this 
paper. The paper, however, will concentrate on the 
evaluation of the adaptive unstructured meshing 
technique for providing solution accuracy as well as 
Increasing computational efficiency for thermal stress 
analysis of structures. 

The governing differential equations for predicting 
the structural response due to both thermal and 
mechanical loads will be presented first. The triangular 
membrane element with temperature distribution over the 
element and a new triangular plate bending element with 
temperature gradient through the element thickness will 
be described. The corresponding finite element 
equations and the associated element matrices will be 
derived and presented. The basic concepts of the 
adaptive unstructured meshing technique and the 
selection of the meshing parameters used for 
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construction of new meshes will be explained. Before 
applying the adaptive unstructured meshing technique to 
complex problems, the technique will be evaluated for 
the stress analysis of a panel with a circular cutout where 
the exact solution Is available. The performance of the 
new triangular plate bending element under the thermal 
load will then be evaluated using an example of a plate 
with temperature gradient through its thickness for which 
the exact solution is also available. Finally the adaptive 
unstructured meshing technique for built-up structures 
will be evaluated by the thermal stress analysis of 
convectively-cooled intersecting frame stiffened panels in 
which the adaptive mesh solution will be compared with 
a fine mesh solution. 


where D » Et 3 /12(1-v 2 ) is the bending rigidity, p(x,y) is 
the transverse pressure, and Mr is the thermal moment 
defined by, 

r 

Mr - Ea (T(z)-T 0 )zdz (6) 
J- V2 


Derivation of Finite Element Equations 

The Constant Strain Triangle (CST) and the 
Discrete Kirchoff Triangle (DKT) finite elements are used 
for the In-plane deformation and the transverse 
deflection, respectively. 


Thermal Stress Analysis Procedure 
Governing Equations 

The equations for the In-plane deformation and 
the transverse deflection of a plate that lies In a local x-y 
coordinate system are briefly described herein. 

In-Plane Deformation. The equations for the in- 
plane deformation are given by the two-dimensional 
equilibrium equations in the form, 

— (Si) + — (Sa} - 0 (1) 

3x dy 


Constant Strain Triangle (CST). The three-node 
CST element assumes a linear displacement 
distribution^ over the element. The element equations 
can be derived by applying the method of weighted 
residuals to the governing differential equation, Eq. (1), 
which leads to the element equations in the form, 

[Kn.HS*,} - (F) + { F m ) (7) 

where the vector {5 m } contains the element nodal 
unknowns of the in-plane displacements in the element 
local x-y coordinate directions. There are two In-plane 
displacements per node or six In-plane displacements 
per element. The element stiffness matrix, [K m ], that 
appears In Eq. (7) is defined by, 


where the vectors (Si) and {S 2 } consist of the stress 
components given by, 

(Si) - [ °x T *y ) T (2) 

(Sd - bxy <VJ T 

The stress components are related to the strain 
components and the temperature by the generalized 
Hooke's law, 

{ 0 } - (CnO{e) - [C,„l{a}(T(x,y)-To) (3) 


where {a} contains the stress components o„, o y , and 
x X y, (C m ] Is the material stiffness matrix, and (a) is the 
vector of the coefficients of thermal expansion. For the 
plane stress case, these material matrices are given In 
Ref. 3. The vector of the strain components Is related to 
the displacement gradients given by, 
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dx dy dy dx 


[Km] - [Bm] T [Cm][B m ]tA (8) 

where the strain-displacement interpolation matrix, [B m ], 
is given in Ref. 3. The first vector, {F}, on the right-hand- 
side of Eq. (7) contains the applied mechanical forces at 
element nodes. The second vector, {F m }, consists of the 
equivalent nodal forces due to the thermal load. This 
second vector is defined by, 

(Fm) - [Bm] T [Cm]{aHWT 0 )tA (9) 


where T«g « J T(x,y)dA/A. 

Discrete Kirchoff Triangle (DKT). The three-node 
DKT element assumes a cubic distribution of the 
transverse deflection^ over the element. The element 
equations can be derived by applying the method of 
weighted residuals to the plate bending equations, Eq. 
(5), which leads to the element equations in the form, 


Transverse Deflection. The equation for the 
transverse deflection in the z-direction normal to the x-y 
plane of the plate is given by the equilibrium equation 4 
in the form, 

_ 4 -.4 -4 

_.dw „ d W d W . 

D ( +2 + — -) - 

dx * dx*dy 2 dy 4 

< 5 > 


[Kb]{6b) = (F p ) + {F b } (10) 

where the vector {5b} contains the element nodal 
unknowns of the transverse deflections and the rotations 
(slopes). Each node has a transverse deflection in the 
element local z-coordinate direction and two rotations 
about the element local x-y coordinate directions. Thus 
there are nine degrees of freedom per element. The 
performance of the element has been evaluated 
thoroughly in Ref. 6 using several benchmark problems 
for bending of plates under mechanical loading. 
However, the performance of the element under thermal 
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loading has not been evaluated. Without thermal 
loading, only the element stiffness matrix, [Kb], and the 
nodal force vector due to the applied pressure, {F p }, that 
appear in the element equations, Eq. (10), are needed. 
These matrices are defined by, 

[Kb] « | [Bb] T [Cb][Bb]dA (11) 


(F p ) = J [N b ] T pdA (12) 


where the strain-displacement Interpolation matrix, [Bb], 
the plate material stiffness matrix, [Cb], and the plate 
element interpolation function matrix, [Nb], are given in 
Refs. 5-6. The above two element matrices, [K b ] and 
(F p ), can be evaluated in closed form (i.e. numerical 
integration is not required) and the details of the 
derivation are given in Refs. 6-7. 

With the presence of thermal loading, the vector 
of the equivalent nodal forces, (Fb), due to the thermal 
moment is included in the element equations, Eq. (10). 
This additional vector is defined by. 


(F b ) - J [Bb] T (M)dA (13) 

where 

(M) « [Mr Mt 0] t (14) 


The thermal moment Mj is defined in Eq. (6). The 
vector of the equivalent nodal forces due to thermal 
loading, {Fb}, above can also be derived in closed form. 
For an element with arbitrary temperature distribution 
through the plate thickness, T * T(z), this vector, {Fb), Is 


{F b } 


Mt[G] t 


(15) 


where the matrix [G] has been derived and given in the 
Appendix. 

Note that, for general built-up structures, elements 
can be arbitrarily oriented in three dimensions. 
Transformation of these element matrices from a local x- 
y-z coordinate system to a global X-Y-Z coordinate 
system is required before assembling them to obtain a 
set of simultaneous equations. In such a global X-Y-Z 
coordinate system, each node has six degrees of 
freedom which are the three translations and the three 
rotations. The matrix transformation procedure, 
however, can be easily worked out and is not presented 
herein for brevity. 


previous mesh. The new mesh will have small elements 
(short nodal spacings) in regions of large changes In 
solution gradients and large elements (large nodal 
spacings) in regions where the gradient changes are 
small. Proper nodal spacings used for constructing a 
new mesh are determined by following a procedure 
analogous to the solid mechanics concept of finding the 
principal stresses, oi and 02 , from a given state of 
stresses, <y x , o y . and t*y, I.e., 



*xy 


Oi 

0 






*xy 



0 

02 . 


At a typical node In the previous mesh, the second 
derivatives of the key parameter for remeshing, <(», 
(analogous to the stress components in Eq. (16)) are 
computed and the two eigenvalues (analogous to the 
principal stresses) are then determined. 


a 2 <t> 

dx 2 

a 2 * 
ax ay 
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The larger eigenvalue, X = max (Xi,X 2 ), Is then selected 
for that node and the same process is repeated for all the 
other nodes in the previous mesh. Proper nodal 
spacings, denoted by h, used for constructing a new 
mesh are then determined from an equidistribution 
condition 1 , 

Xh 2 - constant « Xmuhmin (18) 


where X max Is the largest eigenvalue of all the nodes in 
the previous mesh and h m i n is the specified minimum 
nodal spacing for the new mesh. 

The second derivatives of the key parameter, $, 
are determined using the following procedure. The first 
derivative (e.g. with respect to x) is assumed to vary over 
an element in the form, 


. [ N J |— 1 

dx W 


(19) 


where the vector on the right-hand-side contains the 
unknown nodal values of the first derivatives. These 
unknowns are obtained by solving, 



( 20 ) 


Adaptive Unstructured Remeshlng Procedure 
Remeshlng Strategy 

The basic idea of adaptive unstructured 
remeshing 1 Is to construct a completely new mesh based 
on key parameters from the solution obtained on a 


where a$/3x is the computed uniform element gradient. 
The process is repeated to determine the second 

derivative d^/ax 2 . Determination of the other second 
derivatives which appear in Eq. (17) is performed in the 
same fashion. 
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Remeshlng Parameters 

The adaptive remeshing technique requires a 
selection of proper key parameters (<t> in Eq. (17)). For 
structural problems under mechanical load alone, stress 
is an appropriate choice. However, the key parameter 
representing the stress should be a scalar quantity which 
takes into account all the stress components, such as the 
Von Mises stress defined in two dimensions by, 


OVon Mises “ V (a* - Oy)^ + O? + Oy + 6t?y (21) 

For structural problems under both mechanical and 
thermal loads, the structural temperature and the Von 
Mises stress should be used as the key parameters for 
remeshing simultaneously so that the mesh generated 
can accurately represent the prescribed temperature 
distribution and can capture the high thermal stress as 
well as the mechanical stress concentration. 


Applications 

Three example problems are presented. The first 
example demonstrates the effectiveness of the adaptive 
unstructured remeshing technique on a plane stress 
problem using CST elements. The second example is 
used to evaluate the performance of the triangular DKT 
plate bending element under thermal loading. The third 
example combines both the CST and DKT elements and 
demonstrates the capability of the adaptive unstructured 
remeshing technique for a built-up structure. 


Panel With Circular Cutout 

A panel with a circular cutout, shown in Fig. 1, is 
an ideal structure to demonstrate the capability of the 
adaptive unstructured meshing technique because 
closed-form solutions** are available. The panel is 
subjected to an applied uniform stress (o 0 ) of 10 ksi in the 
longitudinal y-direction (see figure). The exact normal 
stress (oy) distribution along the x-dlrection at y=0 is 
given by, 



where a is the radius of the circular cutout. The peak 
stress (3o 0 ) of 30 ksi occurs at the boundary of the cutout 
(point A In the figure). 

Due to symmetry, a quarter of the panel is 
analyzed. The initial unstructured mesh consisted of 1 84 
nodes and 306 triangles as shown in Fig. 2(a). Details of 
the mesh near the cutout are shown in Fig. 2(b). With this 
mesh, a plane stress analysis with CST elements was 
performed to predict the panel deformation and the stress 
distributions. The contours of the normal stress (o y ) 
distribution In the longitudinal y-direction near the cutout 
are shown in Fig. 2(c). The peak predicted stress (o y ) at 
the cutout is 19.5 ksi compared to the exact peak stress 
of 30 ksi. The 35% error In the peak stress is due to the 
large elements used near the cutout region. 


The Von Mises stresses obtained from this initial 
mesh solution are used as the key parameter for the 
adaptive remeshing. The new adaptive mesh, shown In 
Fig. 3(a), has approximately 20% fewer nodes and 
elements (149 nodes and 250 triangles) than the initial 
mesh. However, small elements are concentrated in the 
region of high stress gradients near the cutout (shown in 
Fig. 3(b)) to provide a more accurate stress solution. The 
normal stress (o y ) contours obtained with this mesh are 
shown in Fig. 3(c). The peak predicted stress (o y ) at the 
cutout is now 28.5 ksi which is only 5% lower than the 
exact level of 30 ksi. 

The normal stress (o y ) distributions along the x- 
direction at y=0 obtained from the Initial and the adaptive 
meshes are compared with the exact stress distribution in 
Fig. 4. By defining a percentage error based on the L -2 
norm of the form, 

( f l(<r y )Exact * (Oy)F.E.) 2 dx) 1/Z 

% error . — x 100 (23) 

(f l(o y )Ex»ct) z dx) 1/z 


the figure shows that the solution error was reduced from 
17% to 2% with one adaptive remesh. In addition, 
increased solution accuracy was obtained with fewer 
unknowns because small elements are concentrated in 
the region of high stress gradients near the cutout and 
large elements are generated in other regions. 

It is important to note that the adaptive meshing 
technique automatically generates refined elements in 
the regions of high stresses. A priori knowledge of the 
solution to the problem (e.g. high stress regions that 
require refined elements) Is not needed before 
performing the analysis. The technique thus provides an 
advantage over the standard finite element procedure 
especially for more complex problems or larger 
structures (such as the structure which will be presented 
in the last example) where a priori knowledge of the 
solution does not exist. 


Plate Bending Due To Thermal Load 

The unconstrained aluminum plate with a 
temperature gradient through its thickness, shown in Fig. 
5, is a simple example for evaluating the performance of 
the DKT plate bending element. The exact solution 4 for 
the transverse deflection is given by, 

w (x,y) - -lML( x 2 + y2) (24) 

4 E r 


where Mr is the thermal moment defined in Eq. (6). For 
a linear temperature distribution through the plate 
thickness with a temperature Tjop on the top surface and 
T Bottom on the bottom surface, the thermal moment is 
given by, 

Mr - i2 0^ T (TTcp ' TBoBom) (25) 
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Due to symmetry, a quarter of the plate is 
analyzed. The finite element model consists of 32 DKT 
plate bending elements and 25 nodes with 5 nodes 
equally spaced In the x and y directions (See Fig. 5.). 
The top and bottom surface temperatures are 630°R and 
530°R, respectively. A comparison of the finite element 
and the exact transverse deflection solutions along the x- 
direction at y=0 is shown In Fig. 6. The figure shows that 
the DKT plate bending elements essentially provide an 
exact transverse deflection solution. 

As described earlier, all the DKT plate bending 
element matrices are evaluated In closed form. This 
simplifies the programming task, reduces the 
computational effort and eliminates numerical error that 
may occur in generating such element matrices. In 
addition, reference 6 has demonstrated the effectiveness 
of this element for providing solution accuracy for plate 
bending problems under different types of the 
mechanical load. The example presented herein 
demonstrates the performance and the capability of the 
DKT element for providing accurate structural response 
due to the thermal load. The next example will combine 
the CST elements (used in the first example) and the 
DKT elements (used in this example) with the adaptive 
unstructured meshing technique for thermal stress 
analysis of a more complex built-up structure. 


Thermal Stress In Intersecting Panels 

To demonstrate the capability of the adaptive 
unstructured meshing technique for thermal stress 
analysis of more complex structures, a three-dimensional 
structure which represents a scramjet engine inlet is 
considered. The structure experiences high heat transfer 
rates from the Impingement of an oblique shock from the 
vehicle forebody as illustrated in Fig. 7. The localized 
heat transfer rates result in high temperature gradients 
and their attendant thermal stresses in the engine panel. 
A typical 'built-up' engine inlet structure, shown in the 
lower right corner of the figure, may consist of panels and 
stiffeners and may be actively cooled. To perform the 
thermal stress analysis of such a structure, a finite 
element model consisting of two-dimensional plate 
bending and membrane elements for the panels and the 
stiffeners is required. Note that, In constructing such a 
finite element model, the same nodal discretization is 
generally required along the Interface between the 
stiffeners and the panels. Such a discretization 
constraint, combined with the requirement of placing fine 
elements at the high thermal stress regions for an 
accurate thermal stress solution, is tedious if done 
nonadaptively. Furthermore, several modelling iterations 
may be needed before achieving the desired solution 
accuracy. 

For the purpose of demonstrating that the thermal 
stress solution accuracy can be improved by using the 
adaptive unstructured meshing technique, both the 
structural temperature distribution and the structural 
boundary conditions are prescribed as shown in Fig. 8. 
The structural boundary conditions are applied such that 
the intersecting panels can move freely in the global X, Y, 
and Z directions as illustrated in the figure. The 
prescribed temperature distribution simulates the 
structural temperature response from the shock 
impingement illustrated in Fig. 7. The distribution is 
assumed to decay radially as a quadratic function from a 


peak of 700°R (at the centers of the two "hot spots') to 
50°R (in the regions away from the two hot spots) as 
illustrated in Fig. 8. The temperature of 50° R is 
representative of the coolant temperature. 

For most thermal stress problems, high 
compressive stresses normally occur at the high 
temperature and temperature gradient regions. Refined 
finite elements are thus needed in these regions. For 
more complex structures, high thermal stresses may also 
occur In other regions that have lower temperature or 
temperature gradients depending on the complexity of 
the structures and how they are constrained. Thus, in the 
standard finite element analysis procedure, a relatively 
fine finite element model is normally employed first to 
identify such high thermal stress regions. An example of 
a fine finite element model Is shown in Fig. 9(a) which 
consists of 3,168 nodes (19,008 degrees of freedom) and 
6,1 1 4 triangular elements. 

The predicted thermal stress distribution in the Y- 
direction superimposed on the deformed geometry 
(greatly exaggerated), obtained from this fine finite 
element model, is shown In Fig. 9(b). As expected, high 
compressive stresses with a magnitude of 25 ksi occur at 
the two hot spots. However, a tensile stress of 54 ksi 
occurs between these two regions. Such phenomenon 
may not be anticipated before analyzing the problem or 
before constructing the finite element model. This 
example demonstrates a need for an automated analysis 
procedure that can provide an accurate solution without 
a priori knowledge of the solution. 

The solution for this fine finite element model 
requires a computational time of approximately 13 CPU 
hours on a VAX-8550 computer. This computational time 
can be reduced significantly by placing larger elements 
in the regions away from the hot spot regions. As 
mentioned earlier, construction of a finite element model 
with different element sizes while maintaining the same 
nodal discretization along the interface between various 
structural components is difficult. Furthermore, it is also 
difficult to construct a finite element model using all the 
four-node quadrilateral elements to satisfy such 
requirements and to achieve smooth transition of 
element sizes from the fine to the coarse mesh regions. 
These difficulties can be alleviated and the analysis 
computational time can be reduced by the adaptive 
unstructured meshing technique described next. 

Application of the adaptive unstructured meshing 
technique starts by constructing a fairly uniform but 
coarse initial mesh as shown In Fig. 10. This initial mesh 
consists of 547 nodes (3,282 degrees of freedom) and 
994 triangular elements, which is approximately 1/6th the 
size of the previous mesh. The predicted thermal stress 
distribution in the Y-direction superimposed on the 
deformed geometry is shown In Fig. 11. The peak 
compressive stress at the two hot spots is 20 ksi while the 
peak predicted tensile stress between the two hot spots 
Is 33 ksi (approximately 40% lower than the fine mesh 
solution shown In Fig. 9). The computational time 
required by this initial mesh is about 15 minutes on the 
same VAX-8550 computer. 

The temperature and the stress distribution (in 
form of the Von Mises stress given by Eq. (21)) obtained 
from the Initial mesh are used as the meshing 
parameters to automatically construct a new adaptive 
mesh. The new adaptive mesh, shown In Fig. 12, has 
slightly fewer nodes and elements (528 nodes and 989 
elements) than the initial mesh. However, smaller 
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elements are concentrated in the high stress regions on 
both the panels and the stiffeners. Larger elements are 
generated in other regions to reduce the problem size 
and thus the computational time. The predicted thermal 
stress distribution In the Y-direction superimposed on the 
deformed geometry is shown in Fig. 13. The figure 
shows that the adaptive mesh model provides the same 
stress solution level as obtained from the fine mesh 
model shown in Fig. 9 (i.e. the peak predicted tensile 
stress of 54 ksi is between the two hot spots that have 
compressive stresses of 25 ksi). The total CPU time used 
for the entire adaptive analysis procedure is 
approximately 30 minutes, I.e. 15 minutes each for the 
solutions on the Initial and the adaptive mesh and 
approximately a minute for constructing the two meshes. 
This total CPU time therefore is 1/25th the CPU time 
required by the fine mesh model. 

The adaptive mesh model shown in Fig. 12 also 
highlights the ease of constructing a finite element model 
by using triangles. These triangles provide smooth 
transition of the element sizes from refined to coarse 
mesh regions, i.e. in the neighborhood of the hot spots as 
shown in the figure. The smooth element transition 
provides a comparable stress distribution to the fine 
mesh solution as indicated by comparing the stress 
contours obtained from the adaptive mesh model (Fig. 
13) and the fine mesh model (Fig. 9). This last 
application demonstrates that, without a priori knowledge 
of the solution, the adaptive unstructured remeshing 
technique provides the same thermal stress solution 
accuracy as the fine mesh solution but at a reduced 
analysis computational cost. 


Concluding Remarks 

An adaptive unstructured remeshing technique 
for thermal stress analysis of built-up structures has 
been developed. The technique generates a new mesh 
based on the solution obtained from a previous mesh. 
The new mesh consists of small elements In regions of 
large change in the stress gradients and large elements 
In other regions where the gradient changes are small. 
The finite element thermal stress formulation for the 
triangular membrane (CST) and plate bending (DKT) 
elements under both mechanical and thermal loads is 
presented. All the element matrices were derived in 
closed form to simplify the programming task, reduce the 
computational time and eliminate error associated with 
the numerical integration. To construct a new adaptive 
mesh, proper nodal spacings are determined from the 
meshing indicators that are represented by the second 
derivatives of the selected meshing parameters. The 
higher the second derivatives computed, the smaller the 
element sizes generated. The meshi ng pa rameters are 
selected according to the type of the problems. Tor 
example, both the temperature and Von Mises stress are 
used together as the meshing parameters for thermal 
stress problems so that refined elements are generated 
In the regions of high thermal stress and mechanical 
stress concentration. 

The two applications presented in this paper 
demonstrate that the adaptive unstructured meshing 
technique: (1) reduces modelling effort because a priori 
knowledge of the solution Is not required; (2) provides 
Improved solution accuracy by adapting the mesh to the 
physics of the solutions; (3) reduces the problem size by 
automatically generating small elements in the regions 


with high solution gradients and large elements in the 
other regions; and (4) provides smooth solution 
transition from fine to coarse mesh regions with the use 
of triangular elements. 


Appendix 


The vector of the equivalent nodal forces due to 
the thermal load, {Fb}, for the DKT plate bending element 
given in Eq. (15) can be evaluated in closed form. This 
vector includes the matrix [G] given by, 


[G] 


yai [Git] + yi2 [G12I 

-X31 [G21] - x« [G22] 

-X31 [G11] - x« [G12] + yai [G21] + yi2 [G22] 


The coefficients xy and yy, i, j»1, 2, 3, are defined in 
terms of the element nodal coordinates by, 


Xy - X,- X| ; y, - yi-yj 

The row matrices, [Gy], i,j=1,2, that appear in the matrix 
[G] above are given by. 


[G,i ] - 

[ P 5 

-ds 

-rs 


P 4 

d 4 

r 4 


-P 4 -P 5 

d 4 -qs 

r 4 -r 5 ] 

[G 12 ] - 

[P6 

-qs 

-rs 


P 4 +P 6 

q4-qs 

r 4 -r 8 


-P 4 

q4 

U ] 

[G 21 ] - 

1*5 

3 -r 5 

qs 


U 

- 3 +r 4 

-q4 


-U-ts 

r 4 -r 5 

-q4+qs] 

[G22] - 

[te 

3 -r 6 

qe 


U+te 

r 4 -re 

-q4+qs 


■u 

- 3 +r 4 

-<U ] 


The coefficients Pk, qk. h* and tk, k = 4, 5, 6 in the row 
matrices, [Gy], above are defined by, 

Pk - -6 xy/lf ; qk - 3 xy yy / 1§ 

rk - 3y*/§; tk - -6y,/§ 


where l|f = x^ + yj and k = 4, 5, 6 for ij = 23, 31, 12, 
respectively. 


Intersecting Panels Dimensions 

The size of the copper intersecting panels used in 
the thermal stress analysis in the last example is 
assumed to be 6 x 10 x 4 inches in the X, Y, and Z 
directions (see Fig. 8), respectively. 
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G 0 = 10 ksl 
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Fig. 2 Initial finite element mesh for the panel with a 
circular cutout and the predicted stress 
contours. 




(a) Adaptive mesh 



(b) Detailed mesh 
near the cutout 


(c) Predicted stress 
contours 
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Fig. 5 Schematic of an unconstrained plate subjected 
to linear temperature distribution through its 
thickness. 
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Fig. 3 Adaptive finite element mesh for the panel with 
a circular cutout and the predicted stress 
contours. 


Fig. 6 Comparison of the predicted finite element and 
the exact transverse deflections along the x- 
direction at y=0 for the plate with linear 
temperature distribution through its thickness. 



Fig. 4 Comparison of the predicted finite element and Fig. 7 Schematic of high heating and thermal stress 
the exact stress distributions along the x- regions on a scramjet engine inlet represented 

direction at y=0 for the panel with a circular by a structure with intersecting panels, 

cutout. 
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(a) Fine mesh 
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(b) Predicted axial stress contours 

Fig. 9 Standard fine finite element mesh and the 
predicted axial stress contours on deformed 
configuration. 
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Fig. 12 Adaptive unstructured finite element mesh for 
the Intersecting panels. 
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Fig. 10 Initial finite element mesh for the intersecting 
panels. 


Fig 13 Predicted axial stress contours from the 

adaptive finite element mesh on deformed 
intersecting panel configuration. 
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